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Abstract 



Long-distance effects in exclusive rare semileptonic transitions B — > (K, K* ) are analysed 
within a relativistic quark model. The meson transition form factors, describing the meson 
amplitudes of the effective weak Hamiltonian, are calculated within the dispersion formu- 
lation of the quark model as relativistic double spectral representations through the wave 
functions of the initial and final mesons. The dilepton spectra and lepton asymmetries are 
considered within the framework provided by the Standard Model. It is found that, while 
the non-resonant decay rates are very sensitive to quark model parameters, the model de- 
pendence of the predicted dilepton forward-backward and lepton polarization asymmetries 
is remarkably small, providing only an overall ~ 10% uncertainty. 
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The investigation of rare semileptonic decays of the B meson induced by the flavour- 
changing neutral current transitions b — > (s, d) provides an important test of the Standard 
Model (SM) and its possible extensions. Rare decays are forbidden at tree level and occur 
at the lowest order only through one-loop diagrams. This fact opens the possibility to probe 
at comparatively low energies the structure of the theory at large mass scales, thanks to the 
contributions of virtual particles in the loops. However, in order to reliably separate (pertur- 
bative) short- distance effects, the (non-perturbative) long-distance contributions, entering the 
amplitudes of exclusive rare decays, should be known with enough accuracy. The theoretical 
investigation of these contributions encounters the problem of describing the hadron structure, 
which yields the main uncertainty in the predictions of exclusive rare decays. In this letter 
exclusive rare semileptonic transitions B — » (K, K*) are analysed within the framework of the 
SM and adopting a relativistic constituent quark (CQ) model. It is shown that the differential 
decay rates are very sensitive to quark model parameters, whereas the model dependence of 
the predicted dilepton forward-backward and lepton polarization asymmetries is quite small, 
providing only an overall ~ 10% uncertainty. 

1. Operator basis. The effective weak Hamiltonian, which describes the b — > s£ + £~ 
transition, has the following form JlJ 

AC 

TUff = -^v*K Oifr) (i) 

where Gp is the universal Fermi constant, the quantities Cj(/i) are the Wilson coefficients, 
obtained after integrating out the heavy particles, and the Oj's are the basis operators; within 
the SM, the operators providing the main contribution to rare decays are 0, |3| 

Ox = - (s Q 7 M (l - 7 5 )6 Q ) (c/37 M (l - 75)0/9) , 

2 = 7 (sa7"(l - 75)&/3) (c/37m(! - 7s)c a ) , 

°7 = ^T^afvM 1 + 7s) + m s (l - 7 5 )]6 Q F^, 
62,71 

O9 = ^(s a ^(l- l5 )b a ) (£^£), 

Oio = ^(s a Y(l-i 5 )b Q ) (£ W £), (2) 

In Eq. (|I]) the renormalization scale \x is usually chosen to be /x ~ in order to avoid large 
logarithms in the matrix elements of the operators 0«. The Wilson coefficients Cj reflect the 
specific features of the theory at large mass scales; they are calculated at the scale /1 ~ M w 
and then evolved down to \x = mb by the renormalization group equations. The analytic 
expressions for Cj(/i) in the SM can be found, e.g., in @. In what follows, the values of the 
Wilson coefficients at the scale \i = = 5 GeV are 0: Ci(m&) = —0.235, C^fm?,) = 1.1, 
C 7 {m b ) = -0.333, C 9 (m b ) = 4.09 and C w {m b ) = -4.32. 
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The four-quark operators 0\ and 2 generate both short- and long-distance contributions 
to the effective weak Hamiltonian (|T]). Both contributions can be taken into account by replacing 



Cgirrib) with an effective coefficient Cg (m b ,q 2 ) given by 



r eff 



(m b ,q 2 



C 9 (m b ) + [3Ci(m b ) + C 2 (m 6 )] 
2 3 



,.„ c q. 
h{ — , — 
nib f^t 



+ 



v, 



E 



nT ( y. _ U )M Vi 
q 2 -iM Vi T 



Ml 



(3) 



where q 2 is the invariant mass squared of the lepton pair. The short-distance contributions are 
contained in the function h(m c /nib, q 2 /ml), which describes the one-loop matrix element of the 
four-quark operators 0\ and 2 (see, e.g., @ for its explicit expression). The long-distance 
contribution, related to the formation of intermediate cc bound states, is usually estimated by 
combining the factorization hypothesis and the Vector Meson Dominance (VMD) assumption 
|3]. !| ; phenomenological analyses [[| suggest that in order to reproduce correctly the branching 
ratio BR(B -> J/^X -> g+t~X) = BR(B -> J/ipX) ■ BR(J/^ -> the fudge factor «, 

which appears in Eq. ([3D to correct phenomenologically for inadequacies of the factorization + 
VMD framework, should satisfy the approximate relation: k [3Ci(m&) + ^(m^)] ~ 1. To sum 
up, the effective weak Hamiltonian has the following structure (cf. [§, §, [|) 



W e// (6 
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where x t = (m t /Mw) an d -^(^t) is given in ||. At m t = 176 GeV one has X(x t ) = 2.02. 

2. Meson form factors. The long-distance contribution to B — > (K,K*) decays is 
contained in the meson matrix elements of the bilinear quark currents appearing in TC e ff, i.e. in 
the relativistic invariant transition form factors of the vector, axial-vector and tensor currents§. 
The amplitudes of meson decays are induced by the quark transition b — > s through the vector 
= 57^6, axial-vector = s^^b and tensor T^ v = sa^b currents, with the following 
covariant structure || 
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b In rare semileptonic decays there is another long-distance effect, known as the weak annihilation, which is 
caused by the Cabibbo-suppressed part of the four-fermion operators not included in the operator basis (0). 
However, the impact of this process in B — > (K, K*) transitions is negligible ||. 
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where q = p\ — p 2 and P = pi + p 2 - 

The relativistic invariant form factors, appearing in Eq. (^|), contain information on the 
non-perturbative aspects of the decay processes, so that they should be calculated within a non- 
perturbative approach for any particular initial and final mesons. To this end various theoretical 
methods have been adopted, like: the light-cone quark model (LCQM) |7|], the constituent quark 
picture f|, the heavy-quark symmetry (HQS) relations |§, the three-point sum rules (3pSR) 



10] and the light-cone sum rules (LCSR) The main outcome of existing analyses is the 

remarkable model-dependence of the predicted form factors, illustrated in Table 1 in terms 
of another set of frequently used form factors, namely: Fi(q 2 ) = f + (q 2 ), F (q 2 ) = f+(q 2 ) + 
q 2 f-{q 2 )/{Pq)i F T {q 2 ) = -(M 1 + M 2 )s{q 2 ) , V{q 2 ) = (M 1 + M 2 )g(q 2 ), A 1 (q 2 ) = f (q 2 ) / \M X + M 2 ) , 
A 2 (q 2 ) = -(M 1+ M 2 )a + (q 2 ) } A (q 2 ) = [q 2 a^q 2 ) + f(q 2 ) + (Pq)a + (q 2 )]/2M 2 , Ti(q 2 ) = -g + (q 2 ), 
T 2 (q 2 ) = -g + (q 2 ) - q 2 g-{q 2 )/{Pq) and T 3 (q 2 ) = (M 1 + M 2 ) 2 [g_(q 2 ) / \Pq) - % 2 )/2]. This fact 
causes quite uncertain predictions for branching ratios and, in particular, for dilepton spectra 
and asymmetries within the SM 0, 0, |9], [10], [11], [1^] . It is clear that such an uncertainty may 



become an obstacle for extracting information on the Wilson coefficients (particularly, their 
signs) and for analysing possible deviations from SM predictions, like those expected in SUSY 
models |, 0. 

We have investigated the relevant meson form factors within a relativistic CQ model 
adopting a dispersion formulation, which has proved to be successful in describing semileptonic 
decays of heavy mesons ||14|| . The dispersion formulation of the quark model has several contact 
points with the LCQM (see, e.g., []15|). However, the LCQM has the problem of a direct 



application at time-like values of q 2 because of the contribution arising from the so-called non- 
partonic diagram, which cannot be killed at q 2 > by an appropriate choice of the reference 
frame. The dispersion formulation overcomes this difficulty. Indeed, the LCQM form factors at 
q 2 < are re-written as double spectral representations in the invariant masses of the initial and 
final qq pairs, and, then, an analytical continuation is performed to reach the time-like region 
q 2 > 0. 

Let us consider the transition from the initial meson q(m 2 )q(m 3 ) with mass M\ to the 
final meson q{mi)q{m 3 ) with mass M 2 , induced by the quark transition m 2 —> mi through the 
current q(m 1 )J^q(m 2 ). For the transition B u — > (K,K*) one has m 2 = mf,, mi = m s and 
m 3 — m u- The CQ structure of the initial and final mesons is described by the vertices Ti 
and T 2 , respectively. The initial 5-meson vertex has the spinorial structure Ti = ^75 G x /\/W c , 
where N c is the number of colours; the final meson vertex has the structure T 2 = 275 G 2 j \[W C 
for a pseudoscalar state and Y 2fl = [A^^ + B{ki — A^)^] G 2 /^/W c , with A = — 1 and B = 
l/(v^2 + Wi + 777,3) for an S*-wave vector meson. At q 2 < the form factors are given by the 
following spectral representation ]TJ| 



00 s t( S2 '9 2 ) ~ 

f( 2 ] _ f ds 2 G 2 (s 2 ) r dsid(si) fj(si,s 2 ,g 2 ) 

m) ~ J n(s 2 -M 2 ) J n( Sl -M?)16\V 2 ( Sl ,s 2 , q 2 ) W 

( m l+ m 3) 2 s-(s 2 ,q 2 ) 

where sf(s 2 , q 2 ) = [^2(777^ + m\ — q 2 ) + q 2 im\ + mf) — [m\ — ml)(ml — m%) ± A 1 / 2 (s 2 , mf, mf) 
A 1 / 2 (g 2 , mf, m2)]/2m\ and A(si, s 2 , S3) = (si + s 2 — S3) 2 — Asis 2 is the triangle function. Equation 
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(Q) corresponds only to the contribution of the two-particle singularities in the Feynman graphs. 
For pseudoscalar and vector mesons with mass M, built up of CQ's with masses m q and nig, 
the function G(s) can be written as 



7r \ s 2 — (m 2 , — rnl) 2 s _ m 2 
G{s) = -^ V , 1 9 q) 3/4 w(k 2 ) (7) 
v2 Js- (m q — rrig) 2 s 

where k = A 1//2 (s, m 2 m?)/2y / i and w(k 2 ) is the ground-state S'-wave radial wave function, 
normalized as / °° dkk 2 \w{k 2 )\ 2 = 1. The double spectral densities /j(si, s 2 , g 2 ), appearing in Eq. 
(||), are explicitly given by 

/+ + /- = 4[mim 2 o;i - m 2 m 3 ai + mim 3 (l - ai) - m 2 (l - ai) + a 2 s 2 ], 

/+ — /_ = 4[mim 2 a 2 - mim 3 a 2 + m 2 m 3 (l - a 2 ) - m 2 (l - a 2 ) + «iSi], 

g = —2/1 [mi«2 + m 2 «i + m 3 (l — ai — cc 2 )] — 45/?, 

a + + a_ = — 4A[2m 2 an + 2m 3 («i — an)] + 4S[Ciai + C 3 a;ii], 

a + - a_ = -\A\-m x a 2 - m 2 («i - 2a i2 ) - m 3 (l - a x - a 2 + 2a 12 )} + AB[C 2 a 1 + C 3 «i 2 ], 

/ = f D + {Ml - si + M 2 2 - s 2 )~g, 

s = 2\m x a 2 + m 2 ai + m 3 (l - a x - a 2 )], 

g + + g- = AA[m 3 (m 1 - m 3 ) + ai(m 1 - m 3 )(m 2 - m 3 ) + a 2 s 2 + 2(3} + 8B(m 1 + m 3 )P, 

g+ — g- = 4A[m 3 (m 2 - m 3 ) + a 2 {m x - m 3 )(m 2 - m 3 ) + a x si] + 8B(m 2 - m 3 )(3 } 

h = -8Aau - 8B[-m 3 ai + (m 3 - m 2 )a u + ("m 3 + mi)ai 2 ], (8) 

where fo = ~^A\m\m 2 m 3 + m 2 (s 2 — m 2 — m 3 )/2 + mi(si — m 2 , — m 2 )/2— m 3 (s 3 — m 2 — m 2 ,)/^ + 
2(3{m 2 - m 3 )} + iBC 3 (3, a x = [(si + s 2 - s 3 )(s 2 - m 2 + mf) - 2s 2 (si - mj + m§)]/A(si, s 2 , s 3 ), 
"2 = [(si + s 2 -s 3 )(s 1 -m 2 + m 2 )-2si(s 2 -m^ + m^)]/A(si, s 2 , S3), P = [2ml-a 1 (s 1 -ml + ml)- 
a 2 (s 2 -m 2 + m 2 )]/4, a n = a\ + A(3{s 2 ) / \{s u s 2 , s 3 ), a 12 = aia 2 - 2/?(si + s 2 - s 3 )/X(s 1 , s 2 , s 3 ), 
Ci = s 2 — (mi + m 3 ) 2 , C 2 = si - (m 2 - m 3 ) 2 and C 3 = s 3 - (mi + m 2 ) 2 — C% — C 2 . 

At q 2 < the representation (|6]) with the spectral densities @ for all the form factors 
but f(q 2 ) coincide with the corresponding LCQM expressions (see, e.g., fl5|l). This is due to 
the fact that within the dispersion approach all the form factors but f{q 2 ) are given by a double 
dispersion representations without subtractions, whereas in order to construct the form factor 
f(q 2 ) from its double spectral density a subtraction procedure is necessary. We fix this procedure 
by requiring that, in case of meson transitions induced by a heavy-to- heavy quark transition, all 
the form factors, including f(q 2 ), satisfy the leading-order Isgur-Wise (IW) relations || as well 
as the subleading 0(l/rriQ) relations of the 1/toq expansion fll6| , which are model-independent 
consequences of QCD. This is fulfilled by the form factors given by Eq. (§) with the spectral 
densities fl8|), provided that the functions Gi(si) are localized near the qq threshold with a width 
of the order Aqcd |17| . Moreover, for meson decays induced by a heavy-to-light quark transition 
the dispersion formulation satisfies the leading-order relations between the form factors of the 
vector and tensor currents given in IJ. 
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The analytical continuation to the time-like region q 2 > generates two contributions: 
the first one is the normal contribution, which is just the expression (Q) taken at q 2 > 0, and 
the second one is an additional anomalous contribution, described explicitly in |14[]. The normal 
contribution dominates the form factors at small q 2 and vanishes when q 2 = (rri2 — mi) 2 , while 
the anomalous contribution is negligible at small q 2 and steeply rises as q 2 — > (772,2 — mi) 2 . 

3. Results. We ran calculations adopting two different CQ models, which will be 
referred to as Set 1 and Set 2. In the former the simple Gaussian ansatz of the ISGW2 model 
is used for w(k 2 ) in Eq. (^), whereas in the latter model the radial function w(k 2 ) is 



the variational solution of the effective qq semi-relativistic Hamiltonian of Godfrey and Isgur 
(GI) ||19|| . The values of the CQ masses as well as the average of the internal momentum 
squared, (k 2 ), are reported in Table 2. The main difference between the radial functions w(k 2 ) 
in the two CQ models relies in their behaviour at high values of the internal momentum k 
(see Table 2 and cf. [|2D], |ST|): while the Gaussian ansatz yields a soft wave function, which 
takes into account mainly the effects of the confinement size, the GI wave functions exhibit 
high-momentum components generated by the effective one-gluon-exchange of the GI potential. 
The impact of these components on the predicted universal IW function and on the form factor 
f+(q 2 ) for heavy-to-light transitions has been analyzed in and found not negligible. 

The results of our calculations of the form factors have been fitted in terms of a simple 
q 2 behaviour of the form }\(q 2 ) = /j(0)/[l — Oiq 2 + o"2<? 4 ]; the resulting values of the parameters 
fi(0), 0i and a 2 are reported in Table 3. It can clearly be seen that the form factors are very 
sensitive to the choice of the quark model. In particular, the form factors obtained with the GI 
wave function (Set 2) are systematically larger than those corresponding to the Gaussian-like 
ansatz (Set 1). This feature is related both to the larger content of high-momentum components 
and to the lower values of the CQ masses characterizing the model of Ref. jHJ with respect to 
the ISGW2 model (see Table 2). 

The behaviour of the form factors of interest is known in case of heavy parent and daughter 
quarks inducing the heavy-to-heavy meson transition Mi —> M 2 . The leading-order I/tuq 
relations read as 

Mi + M 2 Mi - M 2 

J + K ' y+v ; 2^/M^M 2 ^ K ' J K ' y K ' 2y/M x M 2 ^ ; 

s(u) = g(u) = a_(w) = -a+(u) = 2 ^ M ^ 



f(u) = ^MiM 2 (l+u)Z(u), h(u) = -/_= 0(l/m Q ), (9) 

where to = (M 2 + Mf — q 2 ) /2MiM 2 and £(uj) is the IW function. One should not expect these 
relations to work with high accuracy for the B — > (K, K*) transitions, because the daughter s- 
quark cannot be considered heavy enough. As a matter of fact, using the IW function calculated 
in [21], we have checked that large violations of Eq. @ occur, especially far from the zero- recoil 
point. However, Eq. @ implies some (approximate) HQS relations among the form factors, 
namely: Fi(q 2 ) ~ V(q 2 ) ~ A (q 2 ) ~ A 2 (q 2 ) ~ -F T (q 2 ) ~ Ti(q 2 ) ~ T 3 (q 2 ) and F (q 2 ) ~ 
Ai(q 2 ) ~ T 2 (q 2 ), which turn out to be fulfilled within ~ 20% for both Set 1 and Set 2 in the 
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whole kinematical range. 

The averages of our predictions at q 2 = obtained within Set 1 and Set 2 is reported 
in Table 1. It can be seen that our form factors are consistent with the LCSR results. At the 
same time, a striking disagreement with some of the 3pSR results has been found, especially 



in case of the form factors Ft and T3. Indeed, the form factor T3 of Ref. |L0[ has an opposite 
sign and is large in absolute value compared with T 1; as well as the form factor \Ft\ is small 
compared with Fi, at variance with the approximate HQS relations. We want to mention that 
in Ref. |10| both Ft and T 3 are claimed to be constructed in terms of some of the other form 
factors applying equations of motion. However, it can be checked that the expression given in 



10 1 for T3 does not satisfy the HQS requirements in the heavy-quark limit. 

The differential decay rates and asymmetries have been calculated using the expressions 
given in Refs. [10, 12, T3|. The predictions for the dilepton distribution in B — > (K,K*)fi + fi~ 



decays are reported in Fig. 1, where the non- resonant contributions are shown separately. The 
total decay rates turn out to be at least one order of magnitude larger than the non-resonant 
decay rates. However, the resonant contributions are strongly peaked in narrow regions around 
their masses, so that outside these regions the resonance influence is almost negligible. This 
fact allows to reliably separate the resonant contributions from the non-resonant one, which 
contains the relevant information on the Wilson coefficients. In Table 4 our predictions for 
the non-resonant decay rates and branching ratios are listed. The dependence on the chosen 
quark model is strong, yielding a large uncertainty in the predictions; such a drawback may be 
limited by testing the same CQ model in several semileptonic decays, including in particular 
the shape of differential decay rates and lepton spectra, which are expected to be sensitive to 
the form of the CQ model wave functions. Finally, note that the transitions B — > K*fi + ^~ and 
B — > K*e + e~ have different rates, because the amplitude B — > K*£ + £~ has a kinematical pole 
at q 2 = 0, which makes the corresponding decay rate very sensitive to the lower boundary of 
the phase space volume (q 2 = 4m^), while the amplitude B — > K£ + £~ is regular at q 2 = and, 
therefore, insensitive to the mass of the light lepton. 

Our results for total decay rates and branching ratios are summarized in Table 5 and 
compared with the predictions of other approaches. It can be seen that our results are consistent 
with those of Ref. 0, which are based on the application of the HQS relations for the form 
factors, whereas they differ from the predictions of Ref. [JTO] . On one hand, the present level of 
model dependence does not provide the opportunity to extract precise values of V ts and to study 
possible effects beyond the SM. On the other hand, the situation with lepton asymmetries looks 
much more optimistic. Our results for the dilepton forward-backward and lepton polarization 
asymmetries are presented in Figs. 2 and 3. The asymmetries calculated within different quark 
models turn out to differ only by a very small amount (~ 10%) in the whole kinematical range. 
This is due to the following facts: i) the asymmetries depend on ratios among form factors; 
ii) our form factors obey approximate HQS relations, which means that their ratios are only 
slightly model-dependent. We point out that our dilepton asymmetry is approximately a factor 
of 2 lower than the result of Ref. and this disagreement can be traced back again to the 

large difference in the form factor T 3 calculated in |1(| and in the present work (see Table 1). 



c We have taken the resonance phase in accordance with Ref. j|, whereas an opposite sign is used in ]io| . 
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In conclusion, we have analysed rare semileptonic transitions B — > (K, K*) within a 
relativistic constituent quark model, formulated in a dispersion form. The differential decay 
rates are found to be strongly sensitive to the choice of the particular quark model, obtaining 
an overall ~ 50% uncertainty in the predictions. However, at the same time, the asymmetries of 
lepton distributions result to be almost insensitive to quark model parameters, so that they can 
be predicted within the framework of the Standard Model with little uncertainty (~ 10%). Our 
predictions for lepton asymmetries in rare semileptonic decays may provide a reliable starting 
point for investigating possible violations of Standard Model predictions, like those expected in 
SUSY models. 
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Table 1. Transition form factors for the semileptonic decays B — > (K,K*) evaluated at q 2 = within various 
theoretical methods. The last row represents the average of the results obtained in the present work adopting 
the two CQ models, Set 1 (ISGW2 wave function [(is)) and Set 2 (GI wave function jl9||), described in the text. 



Ref. 


fi(0) 


F T (0) 


V(0) 


Mo) 


A 2 (0) 


71(0) 


r 8 (o) 


LCQM p 
3pSR fiqi 
iCS'i? |J 


0.30 

0.25 ±0.03 
0.29 


-0.30 

-0.14 ±0.03 
-0.31 


0.35 

0.47 ±0.03 
0.45 


0.26 

0.37 ±0.03 
0.36 


0.24 

0.40 ±0.03 
0.40 


0.31 

0.38 ±0.06 
0.36 


0.32 

-1.96 

0.39 


T/us work 


0.43 ±0.07 


-0.40 ±0.06 


0.40 ±0.10 


0.34 ±0.08 


0.27 ±0.05 


0.37 ±0.09 


0.35 ±0.09 



Table 2. Values of the CQ masses and of the average internal momentum squared (fc 2 ), in (GeV/c) 2 , for the 
two CQ models adopted in this work. 



Ref. 


m u 


m s 


m c 


m b 


< fc 2 >K 


< fc 2 >K> 


< k2 >b u 


Set 1 


18 


1 


0.33 


0.55 


1.82 


5.2 


0.29 


0.16 


0.28 


Set 2 


20 


] 


0.22 


0.42 


1.65 


5.0 


0.77 


0.34 


0.60 



Table 3. Parameters of the fit fi{q 2 ) = /i(0)/[l — o\q 2 + o^? 4 ] to the B — > (K,K*) transition form factors 
calculated within the two CQ models adopted. 



Decays 


B —f K 


B —> K * 




/+(0) 


/-(0) 


S (0) 


9(0) 


/(0) 


o+(0) 


o_(0) 


MO) 


9+(0) 


9-(0) 


Ref. 














ci 


Cl 


Cl 


Cl 




CT 2 


^2 


CT 2 


CT 2 


CT 2 


c 2 


c 2 


c 2 


c 2 


c 2 


Set 1 


0.36 


— 0.30 


0.06 


0.048 


1.60 


-0.036 


0.041 


0.0037 


-0.28 


0.24 




0.048 


0.050 


0.049 


0.057 


0.0288 


0.053 


0.055 


0.075 


0.058 


0.059 




0.00063 


0.00061 


0.00064 


0.00085 


0.00028 


0.00082 


0.00088 


0.0016 


0.0009 


0.00096 


Set 2 


0.50 


-0.42 


0.080 


0.083 


2.62 


-0.052 


0.067 


0.0085 


-0.47 


0.43 




0.035 


0.036 


0.035 


0.042 


0.011 


0.030 


0.040 


0.044 


0.042 


0.043 




0.00017 


0.00018 


0.00017 


0.00036 


0.00015 


-0.00025 


0.00016 


0.00001 


0.00036 


0.00037 
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Table 4. Non-resonant decay rates and branching ratios of the decays B — > (K,K*) (£ + £ ,uv). Decay rates 
are given in units \V t * s Vtb\ 2 x 10 8 s _1 , while branching ratios are in units \V t * s Vtb\ 2 x 10 -4 . 



Ref. 


B -> 




B -► K*e+e~ 


B -> K*ii + n- 


B —f K*Y, ViVi 


Set 1 


r 


2.0 


25 


5.2 


4.1 


51 






2.9 


37 


7.6 


6.0 


75 


Set 2 


r 


3.3 


39 


11.5 


8.5 


108 




BTZ 


4.8 


57 


17.0 


12.5 


160 



Table 5. Non-resonant decay rates and branching ratios of rare radiative and semileptonic decays of i?-meson. 
Theoretical predictions are given in units | Vt s /0.033| 2 . The results of Ref . Q have been recalculated replacing 
the value | V* s /0.04| 2 with |Vt s /0.033| 2 . The uncertainties in Vt s are not included in the error bars. 



Ref. 


This work 


[3| 


[10| 


Exp. 


Decay 
mode 


r (O 


r is-') 

BTZ 


r (O 

BTZ 


r 

BTZ 


B -> K*~f 


(2.9 ± 1.3) x 10 7 
(3.9 ± 1.7) x 10" 5 


(4.9 ±2.0) x 10" 5 




(4.2 ± 1.0) x 10" 5 p2] 


B -> Kt + £~ 


(2.9 ±0.7) x 10 b 
(4.2 ±0.9) x 10" 7 


(4.0 ± 1.5) x 10" 7 


2 x 10" 7 


< 0.9 x 10~ 5 H 


B — > KY^ViVi 


(3.5 ±0.8) x 10 b 
(5.2 ± 1.1) x 10" 6 


(3.2 ± 1.6) x 10" 6 


(1.6 ±0.4) x 10" 6 




B -> i<T*e + e- 


(9.2 ±3.5) x 10 b 
(1.4 ±0.5) x 10~ 6 


(2.3 ±0.9) x 10" 6 


0.7 x 10~ 6 


< 1.6 x 10- 5 || 


S -> K*fi + n~~ 


(6.9 ±2.4) x 10 5 
(1.0 ±0.4) x 10" 6 


(1.5 ±0.6) x 10" 6 


0.7 x 10~ 6 


< 2.5 x 10~ 5 p4| 


B^K*J2 


(8.7 ±3.1) x 10 b 
(1.3 ±0.5) x 10~ 5 


(1.1 ±0.55) x 10" 5 


(3.5 ±0.5) x 10" 6 
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Figure 1. Differential decay rate for the transitions B — > K£ + £~ (a) and B — > K*/i + n~ (b) versus the invariant 
mass squared of the lepton pair, g 2 , divided by the B-meson mass squared, Mg. Dashed and solid lines correspond 
to our calculations obtained using the two CQ models Set 1 (ISGW2 wave function JOl) and Set 2 (GI wave 
function [p!9|), respectively. The thin lines represent the non-resonant contribution only, whereas the thick lines 
are the sum of non-resonant and resonant contributions. 
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Figure 2. The dilepton forward-backward asymmetry (Afb) for the decay B — > versus q 2 /Mg. 

Dashed and solid lines are as in Fig. 1. The non-resonant contribution (a) and the sum of non-resonant and 
resonant contributions (b) are separately shown. 
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Figure 3. The same as in Fig. 2, but for the lepton longitudinal polarization asymmetry (Pl) for the decays 
B — > K£ + £~ (a, b) and B — > if (c, d), respectively. 
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